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ABSTRACT: We show that two recent independent proposals for regularizing a chiral 
gauge theory stem from one common trick. If the anomaly free complex representation 
carried by the right handed fermi-fields is r one constructs a vector like theory with flavored 
right handed fermionic matter in r + f but with a mass matrix of the order of the cutoff 
and having an index equal to unity in an infinite dimensional flavor space. We present 
a Pauli-Villars realization of the trick that is likely to work to all orders in perturbation 
theory and a lattice version which is argued to produce the correct continuum leading 
order fermionic contribution to the vacuum polarization tensor and readied for further 
perturbative checks. 



INTRODUCTION 



This paper addresses the problem of regularizing a nonabelian chiral gauge theory 
with right handed fermionic matter in an anomaly free complex representation of the 
gauge group. For defmiteness we work in four dimensions and assume the model to be 
asymptotically free so that according to standard lore a continuum interacting theory 
with no free parameters exists. Following [1] we refer to this continuum theory as the 
target theory; a concrete example for a target theory is SO (10) with righthanded fermions 
in the spinorial representation. We shall focus on a particular trick that opens up the 
possibility for a gauge invariant regularization at the expense of introducing an infinite 
set of fermionic regulator fields. Two recent independent proposals [2,3], although quite 
different in appearance, can be built around this trick. 

We start in Minkowski space. The gauge group index will be suppressed and plays no 
role in most of what follows. To each righthanded fermion of the target theory we attach 
a tower consisting of righthanded and lefthanded fermionic fields labeled by a new index 
that lives in a new internal space. The total fermion number of each tower is conserved. 

We wish to give Dirac masses of the order of the cutoff to all the extra fields and, at the 
same time, make the theory look vector like and gauge invariant. We write the Lagrangian 
density as: 

C =Ci + C g (A), Ci = - %gA)^ + ij(MP R + M^P L )^ 

PR=A^, P L= l ~^p L i A li = A°T a , tr(T a T b )=5 ab ^ 

To make the theory look vector like we need to identify the internal spaces associated with 
the lefthanded and righthanded fields and M becomes an operator acting in this space. 
The free fermionic theory will have massless righthanded particles equal in number to tir = 
dim(Ker(M)) and massless lefthanded particles equal in number to ni = dim(Ker(M^)). 
If the internal space has a finite dimension, n^> = n^, and we cannot get to the target 
theory. To make = 1 and = we need an infinite dimensional space; a simple 
example is M = Aa where A is the cutoff and a is a bosonic annihilation operator. 

If M were a finite matrix and M' its restriction to the subspace orthogonal to Ker(M) 
the parity breaking appearance of eq. (1) could be eliminated by tp — > exp^cr^)^ where 
a is a hermitian matrix satisfying exp(2ict) = (M 1 ' M'^)^/^M' '. This doesn't work in our 
case and parity is not conserved. However, one can define U = M^MM^)' 1 / 2 and if 
one "forgets" that U is only formally unitary* the theory looks vector like and easy to 
regularize. 

* \j\\J = 1 but Utf = M^(MM^y l M = 1 - Q where Q is the projector on the zero 
eigenspace of M. 
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One can easily find more examples of M but one should try to keep the structure 
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as simple as possible. Before we had M/A = exp( — 4j-)<9 s exp(^-) where s G R is the 
internal space "coordinate" of the harmonic oscillator. This can be generalized to M/A = 
exp(— F(s))d s exp(F(s)) where exp(— F(s)) is square integrable and the internal space is 
the space of square integrable functions on the line. This choice assures that the kernels of 
M and its adjoint will have the correct dimensions. One can also generalize to an internal 
space corresponding to a complex coordinate with M/A = exp(— F(z, z))d z exp(F(^, z)). 
In the last example M represents one chirality component of a Dirac particle living in 
two Euclidean dimensions in an external magnetic field and the flux of that field must 
be appropriately adjusted to get the correct number of zero modes. The other chirality 
component is represented by M^. 

In all the examples above the internal space corresponds to additional internal dimen- 
sions and ultraviolet behavior might be worse than in the bare theory because there are 
more single particle states at high energies. Therefore, regularization has to deal with 
all the original divergences and also potential new ones. If the spectrum of MM^ is dis- 
crete the fields can be labeled by 0, 1,2, ... and it is natural to view the extra fields as 
Pauli-Villars regulator fields whose masses and statistics are adjusted so that everything 
is ultraviolet finite to any order in perturbation theory. If the spectrum of MM' has an 
unbounded continuous component it is more natural to first try to get rid of the extra 
ultraviolet divergences induced by the internal space and this is easiest done by latticizing 
the variable s and thereby bounding the spectrum of MM^ from above. It is natural 
then to also latticize real space and try to get a non-perturbative definition of the target 
theory with hopes of addressing such long standing vexing basic questions about chiral 
gauge theories as, for example, what their spectrum is. We explore both possibilities in 
this paper. 

We first analyze the Pauli-Villars case and show that it essentially corresponds to the 
suggestion of [2]. Our point of view seems quite different from the one adopted there, 
to the extent we understand it; our final formulae can be made to differ only in some of 
the details and we think that the differences are insignificant. The indications are that 
the method works for an anomaly free situation but breaks down when the gauge group 
is represented by fermions in an anomalous representation. Ungauged, global nonabelian 
symmetries can be preserved only if they are non-anomalous. Next we discuss the lattice 
case and show that it is closely related to the proposal of [3]. Our version will differ from 
this proposal in the manner in which the fermions are coupled to the gauge fields; we 
shall not introduce extra gauge or other bosonic fields beyond those present in the target 
theory. Our study will be somewhat preliminary: we shall provide a necessary tool that 
was missing until now which is an explicit formula for the fermionic propagator and shall 
argue that the vacuum polarization induced by fermions has the right continuum limit. 
We hope to come back in the future with a more detailed perturbative study of the lattice 
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version; as far as we can see now the outlook for this approach is better than anything 
we have seen in years, but a completely well defined construction is still lacking and more 
work is needed. 



PAULI-VILLARS 

We assume that MM^ and M^M have discrete nondegenerate spectra. In addition we 
assume the existence of a "parity" operator S which anticommutes with M, is hermitian 
and squares to unity. The pure gauge part is regularized in some way, unspecified except 
that it preserves ordinary, continuous, space-time structure and gauge invariance. For 
example, a higher derivative regularization, as defined in [4] , does this. We are considering 
in this subsection only perturbation theory and any loop with at least a gauge or ghost 
internal line is assumed to be made finite by the extra regularization. Thus, all we need 
to do is to construct a gauge invariant regularization of the fermionic determinant in an 
arbitrary external gauge field. Not surprisingly, the overall regularization scheme will 
only work if the representation of the gauge group carried by the righthanded fermions is 
anomaly free. 

The free fermion propagator is: 

W - MP R - Mtpj-l = V+ Mt )PR _i_ + tf + M^-j-i^ (2) 

We choose the statistics to be well defined for fields that are eigenstates of S75 . In this 
way the Lagrangian only couples fields of the same statistics: the $ — A term is unity in 
internal space and couples only fields of the same handedness while the mass terms only 
couple fields of opposite handedness and opposite S parity. A fermion loop is given by 
tracing over internal space, spinor index, gauge group index and integrating over the loop 
momentum. Statistics is taken into account by inserting at some point one statistics factor 
equal to — S. To get the right target theory one needs to pick the S parity of the zero 
eigenvector of M as unity. 

As an example consider the lowest order contribution to the vacuum polarization, given 
by < jfi,aji>,b > = ^a&n^j,. After a certain amount of algebra one obtains: 

n A*"(p) = / T^Ns(i-f)- l r^^l l=^^]+\ tr \l-^h^ ^ 

J (2tt) 4 l 2 4+- VmTm 4- - VmTm 2 4+ ^4- J 

q± = k ± ^ and Tr denotes a trace over everything while tr sums over spinor indices. 

The 7 matrix trace in the last term is proportional to e^vaflQ+q-i under k — > —k the 
term is odd and hence the momentum integral over it vanishes. This formal argument* 

* The last term in eq. (3) does not converge by power counting. 
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can be easily made rigorous. Once the first term in eq. (3) is arranged to be finite the 
transversality of the vacuum polarization tensor will be assured by the usual (no longer 
formal) argument. 

To understand finiteness let us consider a slightly more general type of loop integral 
/ ^J(fc) with 

n l Q 

I(k) = polynomial( qj )Tr JJ ( 2 _ M ] M ) ^ ~ o") ( 4 ) 

j=l q j 

where qj = Qj-i +Pj, j = 2, 3, ....n q\ = k and the pj are external momenta that sum to 
zero. The trace is over the internal space. Pick the basis of the space as |0 >, \ E > with 
M^M\0 >= and M^M\E >= E\E >, E > 0. We assume that there are no degeneracies. 
Since S commutes with M^M we have S\E >= rj E \E >, S\0 >= |0 > and r] E = ±1. After 
exponentiating each of the denominators with the help of a parameter tj and rotating to 
Euclidean space it is clear that the question of convergence is settled by the behavior of 
the trace when all tj : s are scaled to zero. Let = r ! we have to consider then the 

behavior of 

Sum = 2^r] E exp(-£ , r) + 1 (5) 
E 

Because of gauge invariance and absence of anomalies we only need to tame the logarithmic 
divergence of a few diagrams. It is then sufficient to make the expression in eq. (5) vanish 
linearly when r goes to zero. This is easily achieved with M = Aa and S = (— l) ata leading 
to Sum = tanh(A 2 r). Better control in the ultraviolet is achievable; for example one may 
be interested in the determinant of $ — y4 in higher even dimensions. An arbitrary degree 
of divergence can be tamed by choosing M = AV aat a and the same statistics operator as 

A 2 

above. Now the sum becomes 6>4(0|e T ) and the Jacobi identity ensures that it vanishes 
faster than any power when r approaches zero. This latter choice has been made in [2]. 

At any order we shall have to deal with one extra piece, similar to the one we discarded 
in the calculation of the vacuum polarization, before we can restrict our attention to sums 
of the kind shown in eq. (4). At higher orders these extra terms converge by power 
counting but diagrams with three and four vertices may be divergent and untamable. Our 
investigations indicate that if the theory is anomaly free no problems are encountered and 
everything goes through as for the vacuum polarization graph. 

Essentially, the Pauli-Villars regularization removes all divergences in the parity con- 
serving part of the fermionic determinant but does not touch the parity odd part of the 
logarithm of this determinant. The latter must be finite and unambiguous "by itself" with- 
out any help from the Pauli-Villars regulator fields and this happens only in the anomaly 
free case. 
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It is important to make sure that the counter terms will not affect the zero mode 
structure of the mass matrices; otherwise the target theory would not be attained. This 
is easily checked to leading order and turns just into a matrix generalization of the well 
known dm j oc my rule. Here this works is an obvious way because the internal space 
is totally decoupled from the real space and the spinor and group indices. However, the 
mechanism that protects the desired zero mode structure from radiative corrections is a 
deeper one and this becomes an essential feature when we turn to the lattice where such a 
decoupling cannot be achieved because of the well known "lattice fermion doubling" . The 
lattice version therefore also avoids infinite amount of fermion mass fine tuning [5] . 



LATTICE 



A simple choice of the mass matrix in the lattice case is based on the "wall" realization 
of M. Its continuum form is M/A = d s + f(s), f(s) = F'(s) with / increasing monoton- 
ically from a negative asymptotic value to a positive one. In particular, this means that 
the internal manifold is not compact. M will have a zero eigenfunction and M t will have 
none; this feature is stable under deformations that do not change the large \s\ behavior 
[6]: Milft and M^M have identical spectra for nonvanishing eigenvalues and therefore any 
small deformation of M cannot move the zero mode of M^M up without pairing it up 
with an eigenstate of MM^ of exactly the same energy; however we assume that there is 
a finite gap between the zero mode and the rest of the spectrum and a small deformation 
cannot drag down an eigenstate of MM^ far enough. It is hoped that radiative corrections 
will be well behaved in this sense if the theory is anomaly free. Plugging the "wall" form 
of M into eq. (1) leads to the interpretation of s as a fifth dimension, the derivative part of 
M and (i.e. M^Ml) combining with the kinetic energy into a five dimensional Dirac 
operator, and the f(s) dependent part (i. e. M+Ml) becoming a variable mass term. 

The spectrum of MM^ is unbounded and to cure this we replace inner space by an 
infinite lattice in both directions, s now denotes the discrete integer labeling the sites 
on this internal lattice and M is replaced by a first order difference operator: M s s > = 
$8+1,8' ~ a ( s )$s,s'- ^ a ( s ) approaches a+ with \a+\ < 1 when s — > oo and a- with |a_| > 1 
when s — > — oo the spectral properties we need are preserved. For definiteness we choose 
a°(s) = 1 — mosign(s + \) with < mg < 1 and denote the corresponding M by M°. 



Putting an Euclidean version (with hermitian 7 matrices) of eq. (1) on the lattice we 
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obtain: 

-C E = - (3 mr[U^x)U v {x + fi)ul(x + v)ut{x)]} + ]T ^(x) lt ,[U^x)^(x + fi) 
- Ul(x - (i)il>(x - fi)} + ^i>{x)[P R {M{Vu)^){x) + P L (M\Vu)^)(x)} 

X 

{Vui>){x) = Y^[U»(x)1>(x + + Ul(x - ijlMx - /jl) - 2il>{x)\ 

(6) 

For noninteracting fermion fields in Fourier space, M depends parametrically on the lattice 
momentum p (— n < p^ < n ) with Vjj = \ = V = — 4 sin 2 ^ so that all the unwanted 
potential zeros of the propagator are moved up in energy by Wilson mass terms. We have 
to require that for V in a region close to p^ = dim(Ker(M (V))) = 1 but for V < — 4 
dim(Ker(M(V))) = 0. We also require that have no zero modes in regions around the 
momenta p with sin(p^) = for any \x. We need to change the asymptotic behavior in s so 
that as V varies over the Brillouin zone the spectra of M(V) and M^(V) change as follows: 
for small momenta that satisfy —V < h (an example for h will be provided later), M has a 
normalizable zero mode and the spectrum of M(V)M^(V) is continuous with a finite gap 
separating it from zero. On the boundary of this small momentum region, —V = h, the 
spectra of M(V)M'(V) and M' (V)M(V) are identical, continuous and gapless. Outside 
the region, where -V > h, the spectra of M (V)Aft(V) and Aft(V)M(V) are identical, 
continuous and there is a gap. All unwanted "doubler modes" of the naive fermionic 
action are in this third region. Note that the abrupt disappearance of the zero mode at 
the boundary between the two regions is a necessary consequence of the generic structure 
of the supersymmetric quantum mechanics governing the internal space and not of the 
specific choices we have made. In short, the trick in the lattice version amounts to walling 
a small region of momentum space around the origin off the dangerous areas inhabited 
by doublers. The spectrum has a nonanalyticity in momentum space as a result and one 
has to check carefully for effects caused by this*. Clearly, these nonanalyticities come 
about because internal space is infinite. The divergence responsible for the singularity on 
the four dimensional Fourier torus is an infrared problem from the point of view of the 
extra dimension. Other approaches [7,8] that tried to use singularities fail because gauge 
invariance forces compensating singular changes in the gauge fermion vertices [9,10,11] 
and these introduce unwanted doublers, ghosts, Lorentz violations or uncomputable non- 
perturbative effects [11,12,13] into the theory. Somehow, this new trick has to avoid all of 
these traps. 

We know that the zero modes of M live in the vicinity of the wall in the five dimen- 
sional world and that the physics we are interested in is the defect dependent part of the 

* Of course, some singularity is needed to avoid doublers in a non-dynamical way. 
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action. To eliminate five dimensional bulk effects the effective gauge interaction induced 

g+ (JJ)+S~ (U) 

by fermions, S e g(U,wall), is replaced by S c q(U) = S c ^(U,wall) cff ^ where 

the S^q{U) come from systems with constant mass terms equal to img. The needed 
subtractions can also be interpreted as representing ghost fields. In analytical computa- 
tions, both in weak and in strong coupling expansions, the needed subtractions should be 
implementable order by order. 

We make our remaining choices so that the homogeneous systems which provide the 
subtractions be as symmetric as possible. We would like to have reflection positivity in 
the five dimensional system in all directions in the homogeneous case and five dimensional 
cubic symmetry when U = 1. Clearly, the relevant five dimensional systems have the same 
Lagrangian as above, only the matrix gets replaced by with the corresponding 
aP(s) replaced by a ± (s) = 1 ± uiq. The five dimensional symmetries we want are then 
obtained by making in each case M(Vj/) = M a + with a = 0, ±. It is easy to see that 
with this choice the origin of four dimensional momentum space is correctly walled off with 
< h = 2mQ < 2. Note that we have more or less naturally ended up in the r = 1 case 
of the Wilson mass term (for notations see [14]); this case has better positivity properties 
than cases where \r\ ^ 1 because it has a positive bounded transfer matrix for single step 
translations [15,16,17] so we should feel no need to generalize [3,18,19]. 

The candidates for observables in the target theory are obtained as follows: Pure gauge 
field operators are naively transcribed. Local four dimensional fermions are represented 
by a sum over s of the bare fermion fields with a weight that samples a finite but wide 
region around s = 0. The width of the weight functions should be taken to infinity at the 
very end. Composite operators of the target theory are best constructed directly in the 
regularized version, for example, the currents associated with global symmetries will have a 
sum over all s in their definition via Noether's formula. With this definition the argument 
that correct chiral anomalies are obtained [3,19,20,21,22,23] applies. The contributions 
of the subtraction terms S^(U) to anomalies cancel against each other because of the 
opposite signs of the mass terms. One expects currents associated with anomalous global 
symmetries to be affected by the infinite sum over s in a non-trivial way, enabling the 
charges to flow out to infinity and hence be nonconserved in any finite five dimensional 
slab containing the defect. This shows that is is unclear whether exact unitarity holds 
if one looks only at such arbitrarily large, but still finite, slabs. A proper construction 
will define the whole theory in such a slab with boundary conditions is the s direction 
(the other directions are irrelevant and can be kept infinite or made finite by standard 
periodic/antiperiodic boundary conditions) chosen in such a way that anomalous charges 
will be able to flow out of the system. If the gauged group is not anomalous one would 
like then to see that a finite limit is obtained when the slab's thickness is taken to infinity. 
The choice of good boundary conditions needs care because in a finite dimensional internal 
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space the zero modes will typically appear both in Af and in Aft; this might be avoidable if 
the boundary conditions replace Af and Aft by matrices that are not each other's adjoint 



but this would give up the unitarity of the theory for any finite width of the slab. 

The advantage of the present approach is that one can set up a meaningful test of the 
construction within perturbation theory and this test can be carried out before the issue of 
boundary conditions is settled since the subtractions remove large s divergences and one 
can work always with an infinite internal space. As in the Pauli-Villars case we consider 
the fermion contribution to the vacuum polarization; this only tests the correctness of 
S c q(U) but is known to be a test failed by all previous proposals. From the imaginary 
part of the polarization tensor at small momenta one can read off the spectrum of fermions 
as seen by the gauge bosons and this is the relevant question. It is obvious that naively the 
method works because the propagator has the right structure at nearly zero momenta by 
construction. The only way this argument can fail in perturbation theory is that there are 
further singularities in the integrand of the relevant lattice Feynman diagrams and they 
make contributions to the continuum limit. If we think about the propagator in a basis 
of eigenstates of Aft(V)Af (V) and Af (V)Aft(V) there is a source of trouble as discussed 
before. We aim to show now that if one represents the propagators in four dimensional 
lattice Fourier space and in s-space one can perform the needed subtractions easily and no 
singularities beyond the "good" ones appear. In particular, none are induced by carrying 
out the s-space traces. 

The propagator is given by a formula similar to eq. (2), now on the lattice and in 
Euclidean space. With p^ = sin(p (U ) and p^ = 2 sin ^ we have: 



Gr(p) and Gi(p) are inverses of second order difference operators and an explicit expres- 
sion can be obtained by standard techniques. Both are symmetric under the interchange 
of s and s' and for s > s' they are given by 



H5>MP/« + M\p))P R G R {p) + H J>mPm + M(p))P L G L (p) 




(7) 



[Gr] ss '(p) 



[Ar(p) 
Ar(p)c 

[Ar( P ) 



- B{p)]e- a +(P^ s+s '+ 2 } + B(p)e- a +^ s -^ for s' > 
a + (p)[s+l]+a-(p)[s'+l] for s > 

- C(p)]e- a -^ s+s '+^ + C(p)e- a -^ s - s '} for s < 



-1 



l>s' 



1 



(8) 
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where 



r 1 1 p 1 1 o 

< a±(p) = arccosh^ ~[a±(p) H H --] k a±(p) = 1 =F m + -p 

^2 a+ » a+ » ) 2 



a-(p)e a -iP) -a+(p)e- a +(p)' a+(p)e a +(p) - a_(p) e - a -(p) (9) 



2a+(p) sinha+(p) ' ' ' 2a_(p) smhct-(p) 

The propagator needed for the subtraction terms are obtained by replacing all a±, a± by 
either a+, a+ or a_, a-. Under each of these substitutions Ar(p), A^(p), B(p) and C(p) 
all become equal to each other. It is easy to see that a±(p) > for all p. Subtractions 
make the sum over s, s' converge exponentially and uniformly in the loop momentum k and 
external momentum p. All ordinary manipulations that prove the transversality of Ii^ v {p) 
and its symmetry properties Yl fll/ (p) = U UfM (p) = IT /ijy (— p) hold after the subtractions. All 
vertex factors are singularity free throughout the Brillouin zone. The only singularities in a 
typical integrand then come from the possible poles in the amplitudes A R {p), A L {p), B{p) 
and C(p). Since we have arranged for M(p) to have a zero mode, there is a singularity at 
p = in Al(p). It is easy to see from the expressions for the amplitudes that this is the 
only singularity. A simple calculation gives 



mo(4-m« 

AlKP) = o — + regular terms (10) 

4p z 



The residue of the pole is just the normalization factor associated with the zero mode at 
p = 0. The usual continuum expressions should therefore emerge. We end our perturbative 
analysis remarking that all ingredients for the calculation of Hfj, u (p) are explicitly present 
and a meticulous calculation will make the above argument foolproof. 

We can also address the issue of how the approach presented here avoids the Nielsen- 
Ninomiya [24] theorem. To this end, consider integrating all the fermion degrees of freedom 
except the right handed fermion at s = 0. The kinetic energy term will be non-local and 
can be extracted from the expression for the propagator at s = s' = 0. We find that the 
kinetic energy term is given by 

By eq. (10) A L (p)p 2 goes to the normalization factor of the zero mode and eq. (11) has 
the expected zero at p = 0. The factor p 2 in the denominator generates fifteen poles 
in the Brillouin zone. Normally these poles would have generated ghost contributions 
[10] and destroyed the theory. But, in this approach, integrating out the fermion degrees 
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of freedom also generates pure gauge terms and these cancel the contributions from the 
ghosts. In earlier attempts to solve the problem [8] , it was understood that the bgauge part 
of the action has to be suitably adjusted to cancel the ghost contributions. The approach 
presented here can be thought of as one way to systematically achieve the apparently 
difficult aim faced before. 

FINAL COMMENTS 

It should be clear from this paper that there is a large amount of latitude in choosing 
M and the associate internal space. It is not clear at the moment whether sticking to 
the "wall" situation is the best, but certainly it seems a worthwhile route to pursue, in 
particular because it gives a familiar picture for where the anomalous charge deficit is 
going. The basic trick however may be useful even if the "wall" realization turns out to 
fail ultimately; in its essence what the trick does is to provide the "infinite hotel" that has 
been argued to be a necessity for the existence of genuinely chiral fermions [24] . 
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